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Abstract. We consider social systems in which agents are not only characterized by their states but also 
have the freedom to choose their interaction partners to maximize their utility. We map such systems onto 
an Ising model in which spins are dynamically coupled by links in a dynamical network. In this model there 
are two dynamical quantities which arrange towards a minimum energy state in the canonical framework: 
the spins, Si, and the adjacency matrix elements, Cij. The model is exactly solvable because microcanonical 
partition functions reduce to products of binomial factors as a direct consequence of the Cij minimizing 
energy. We solve the system for finite sizes and for the two possible thermodynamic limits and discuss the 
phase diagrams. 

PACS. 87.23.Ge Dynamics of social systems - 89.75.Fb Structures and organization in complex systems 
- 05.90.-|-m Networks and genealogical trees - 75.10.IIk Classical spin models - 75.10.Nr Spin-glass and 
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1 Introduction 

Properties of many statistical systems are not solely char- 
acterized by the states of their constituents, but also de- 
pend crucially on how these interact with each other, i.e. 
their network (linking) structure. The way networks func- 
tion can often not be fully understood by their linking 
structure alone because function may depend heavily on 
the internal states of individual nodes. Especially social 
and economical interactions are of this kind. Not only ac- 
tions (states) matter but the possibility of choice with 
whom to interact (linking) plays a crucial role in socio- 
economical dynamics [T] . It is therefore tempting to study 
the co-evolution of network structure and internal states. 
In the simplest case, this can be done in the framework 
of the Ising model, which immediately reminds of spin- 
glass models, such as the SK-model [5] or random-bond 
models, see e.g. Ising models where both, spins and 
interactions, are governed by dynamical rules have been 
studied assuming different timescales of evolution, where 
typically interaction topology 'slowly' adapts in a pre- 
determined way on 'fast' relaxing spins [4l. Recently, such 
systems have been analyzed with the replica approach in 
the grand-canonical ensemble assuming that the interac- 
tion topology also minimizes the energy of the system |5] ; 
the coupling of both, spins and interactions to heat-baths 
at different temperatures can be treated in the respective 
formalism as well Note that these doubly-dynamical 
models are in marked contrast to the Ising model on fixed 
network structures, see e.g. [9]. Complementarily the for- 



mation of network structure driven by various Hamilto- 
nians has been investigated in some detail [TU]. We think 
that a full understanding of many processes taking place 
in networks can only be achieved in a combined approach. 
In the following we show that a spin system in the canon- 
ical ensemble where both, linking structure (given by the 
adjacency matrix Cij) and spins Si, minimize the energy, 
can be exactly solved since partition functions reduce to 
products of binomials. 

The following model is classically phrased in terms of 
magnetization of Ising spins. However, the main idea is 
that it can be one-to-one related to economic terminol- 
ogy. Magnetization, m, correspond to market shares in a 
situation of a two-company world. Think for example that 
there exist two telephone providers, A and B. The monthly 
cost for each individual depends on its local connectivity 
(telephone call network) and on the costs per call (intra- 
provider and out of networks calls) fixed by the provider. 
Here the state of an individual, s^, being customer of com- 
pany A would relate to spin up, customers of firm B relate 
to spin down. Connectivity, , is determined by who calls 
whom. It is assumed that fully rational agents minimize 
their costs. The amount of rationality is modeled below 
by temperature, T. The external field, h, in the following 
relates to external biases, such as asymmetries of PR ac- 
tivity in the firms. There are no conceptual problms to 
extend the methodology of the present work to e.g. the 
Potts model, reflecting a more realistic situation of multi- 
ple companies on the market. 
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2 The model 



We study the Hamiltonian 



i>j 



(1) 



where sums are taken over all N nodes of the system. The 
position of the links in the adjacency matrix Cije{0, 1} is 
a dynamical variable. The system has thus two degrees of 
freedom both minimizing energy: the orientation of the in- 
dividual spins Sie{— 1, 1} as usual, and the linking of spins, 
Cij. Cij = 1(0) means nodes i and j are (un) connected. 
We consider undirected networks {cij — Cji), the case of 
directed networks is a trivial extension as pointed out be- 
low. We denote the number of spins pointing upward by 
n| = J2i the number of links hy L = J2i<j '''ij^ mag- 

— ' connectivity c — — 



netization m 



N 



TV' 



and connectedness tp — In the grand-canonical en- 
semble this Hamiltonian was studied in [5j by use of the 
replica method. In this work, we limit our interest to the 
canonical framework. 

We start our analysis with the microcanonical parti- 
tion function for energy E 

N 
N 

= n{N,n^)Q{E,N,L,h,n^) , (2) 

—0 

where f2{N, n|) is the number of configurations for a given 
rij. f2{E, N, L, h, n|) denotes the microcanonical partition 
function for a fixed n| . 

In Eq.Q the calculation becomes greatly simplified 
when realizing that a fixed number of spins pointing up- 
wards, n|, alone is sufficient to determine the spin-state 
of the system since one deals with all the different topolo- 
gies for a given value of . In other words, the crucial 
observation is that the exact spin-configuration {s^} loses 
its relevance because the topology of the network is not 
fixed. In this case partition functions simply reduce to bi- 
nomial factors. 



N 

E 



= 2 



TV 



(3) 



and the remaining task is to determine il{E, N, L,n^). To 
find the number of microstates leading to energy E for 
fixed n|, the only relevant physical fact is whether a link 
£ connects two spins of (un)equal orientation, thus con- 
tributing a unit — J (J) to total energy. The possible en- 
ergy states are Ee{—LJ—Nhm, —LJ+2J~Nhm, LJ— 
2 J — Nhm, LJ — Nhm} where the lowest energy —LJ — 



Nhm is realized if all links connect spins of equal orien- 
tation. In general, if k links connect spins of equal orien- 
tation (L — k links connect spins of different orientation) , 
E — LJ — 2k J — Nhm. It is easy to see that the number 
of possible 'positions' of linking spins of equal orientation, 
Ce, and unequal orientation, a„, is given by 



ae(^,n|)= i(n|(n|-l) + (7V-n|)(iV 
a„(A^, n|) =n|(7V-n|) , 



1)) 



(4) 



for undirected networks. Directed networks trivially follow 
bom af'{N,n^) = 2ae{N,n^) a.ud at" {N , n^) = 2a„(iV,n|), 
because while in the undirected case, < L < N{N—1) /2, 
in the directed case we have, < L'^" < N{N — 1). Each 
link positioned in ae(„)(iV, nt) contributes — J (J) to the 
total energy E. Given Eq. (Hi), the microcanonical parti- 
tion function for given rij and the total partition function 
read 

i^[Jij,l\,L,n,n^) - \ ,^Lj_E-Nhm) ) [ (LJ+E+Nhm) ) ' 

(5) 

U[E,I\,L,h) - I ) ( lLJ^E~Nhm) j [ {LJ+E+Nhm) 



2J 



2J 



(6) 

We can now directly approach the problem of calcu- 
lating the canonical partition function Z(T, N, L, 7i|) of a 
system with fixed rij via the Laplace transform. 



(7) 



Performing the energy summation the exact solution is 

Z(/3, TV, L, /i, n.) = e"f'+'"'"'^— , 

' ' ' r{l + L)r{l + a^- L) 

X 2<^i(-au, -L, 1 + fle - e'^^P) , (8) 

with 2^i(-a,6, -c,a;) = ELo ^'^(-i)^J\' hypergeo- 
metric function and the Gamma function r{x). The total 
canonical partition function finally is 



TV 



Z{(3,N,L,h)= J2 



N 



Z{P,N,L,h,n^) 



(9) 



and all thermodynamic quantities of interest are given ex- 
actly for finite sized systems, of (fixed) dimensions L and 
N . In Fig. [l] we show the internal energy U and mag- 
netization as a function of temperature for different val- 
ues of connectivity c as calculated from Eq. Perfect 
agreement with Monte Carlo simulations of finite sized 
systems is found, where rewiring and spin-fiipping have 
been implemented by the Metropolis algorithm. We note 
that for low connectivities the obtained solutions are in 
very good agreement with the result of independent spins, 
i.e. U = ctanh(/3), as expected. 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 



3 



Z)-0.6 




Fig. 1. Internal energy for = 100, and connectivities c — 
1,3,6. Solid lines correspond to the exact finite size solution, 
Eq. ([9|. Symbols are the results from a Monte Carlo simulation 
of the canonical ensemble pertaining to the Hamiltonian of Eq. 
([T]). Inset: maximum of magnetization, m as a function of T 
(lines exact, symbols MC). 



3 Thermodynamic limits 

Assuming large N, with Stirling's approximation in the 
form (l) - [(&/a)''/° (1 - °> and the notation 

y = £j,F.q. ^ reads 



Z = 2^-^NL{2ifi)-^ J dm J 



dm / dy[I{m,y,L,N)] 



with 



(10) 



/(to, y, L, iV) = exp(-/3Jy) (1 - TO^) 2c 

ml y 

1 — to\ ^ /I — TO^\^ f — — TO^)^ ^ 



1 + TO 



X ( 1 - 

1 + TO^ 



X 1 - 



1 — TO^ 



(l + y)(l + TO2) 

-4t;(i+™'-2¥>(i-2/)) 

-4^(l-'"^-2v(l+2/)) 



(11) 



where c and ip are shorthands for L/N, and L/N"^ , respec- 
tively; h is set to zero for simplicity. In the thermodynamic 
limit Z is reasonably approximated by the maximal con- 
figuration, i.e. the solution to dl /dy — 0, which is 

_ -1 - TO^t + ^(1 -t- TO2t)2 _ 8(p(to2 - {2(p - l)t)t 

(12) 

with t = tanh(/3J). The other solution is outside the al- 
lowed parameter region of y. To ensure a real valued par- 
tition function the conditions, 1 > 2ip , and, 1 > 

2y) x-m? ' have to hold. Regions where they do not hold are 
forbidden zones in the y — m plane, where the integrand 
of Eq. ( 10 1 is not defined, see Fig. [2] It can be shown that 
the maximum condition line, ?/max("T'), always stays in the 
allowed zone, V -1 < to < 1, < t < 1, and < < 1/2. 
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Fig. 2. Logarithm of Eq. ( |11| in the m — y plane for (p — 0.4 
and c = 20000. The forbidden zones are clearly visible. The 
maximum is always reached in the allowed zone. 



For the following discussi on l et us compute the deriva- 
tive of (the log of) / in Eq. (Ill, with respect to m 



2^T 



dlog(/) m ^ 1 

dm " 2^ 1-2^^ 



1 1 — TO 

— log 

2c 1 -t- TO 



1 2^{l-y^)-il-y){l^m^) ' 

2 ^2^(l-y2)_(i + y)(i + ^2) P-J^ 



y',(13) 



ith 



I m ( 1 + m"^ — Aip \ 

^ ~ 2^ \ ^(1 + mH)^ - 8(^(to2 - {2ip - l)t)t ~ J 

(14) 

It is now natural to consider two distinct thermody- 
namical limits, N oo, one, by keeping connectivity c = 
L/N, the other by keeping the connectedness (p = L/N'^, 
fixed. 



3.1 c=const. limit 

We fix c and take N ^ oo. Consequently , ip vanishes as 
1/iV, and the maximum condition from Eq. (12) reduces 
to 



lim J/max 



t + m^ 

1 + TO2t 



(15) 



The limiting cases for infinite and zero temperature can 
be worked out immediately. 



3.1.1 The low temperature case, (3^1 

For Jf3 S> 1, t — > 1, and the maximum condition further 



simplifies to, y„ 



4, and y'^ 



Using this in Eq. (13), and setting ^^^^ 

1 — TO / 4cm 
= exp 



2m f t(t 

- yields 



1 



(16) 
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(b) 
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Fig. 3. Fixed-c thermodynamic hmit. (a) Magnetization m 
as a function of connectivity c at zero temperature for J = 1. 
Below c = 1/2 there is no possibility for magnetization in the 
system. At c ~ 4, practically full magnetization is reached, (b) 
Phase diagram in the T — c plane for J = 1. 



The self-consistent solution is shown in Fig. |3ja): We find 
zero magnetization below a critical connectivity c < 1/2, 
as well as a region where m ^ 0. One can show [TT] that 
Eq. ( 16 1 can be obtained from the results pertaining to 



the grand-canonical ensemble studied in [5] in the limit of 
an infinitely large chemical potential. 



3.1.2 Infinite temperature, (3-^1 

For J/3 <C 1, the maximum condition here becomes, j/max = 
—m?, and y^j^x ~ —^m. Proceeding as , from c?log(/) /dm 



reduces for m ^ to the critical line given by 



0, we get log (^Yq:^ j — —^cm(3J, or m — tanh(2cm/3 J). 

This implies no magnetization for small /3, or T — )■ oo, 
Vc. The phase transition line, separating the phases m = 
and 771 7^ 0, is found by noting in Eq. (Ill that for 
c = const. 

limjv^oo lim, y,L,N)^ e^^-PJy) ("i _ - k 



Differentiating Eq. (17) w.r.t. m, and setting it to zero, 
i.e. 







— log — 
2c ^ 1 



2m 



y 



rn 



1) 



log 



1-y 



log 



1 -m4 

1—771^ 
1 + 771^ 



(1 + trflty 



2/3J 



(18) 



1 



= tanh 



J 

hm,. 



(19) 



-t. The 



Here we used lim.,„^o Vmi 
phase diagram is shown in Fig. |3]^b). 

Let us finish the c =const. case with a statement on 
the critical exponent By inserting Umax, given in Eq. 
(p5|, into Eq. p^, we get 



l+m 



log ft! 

Amc — m? log 
> 0, where T = 



l+m 
1 — m 



(20) 

1//3. Using the Ansatz 



Define r 

777, = Xt^ , where A is a free parameter, in Eq. (20 1, for 
T <C 1 one gets the expression 

120(1 



4c2 , 



1 



J_ 

2c 



(21) 



Since the left hand side does not depend on t, the only 
possible choice for the critical exponent is (3 — 1/2. This 
is exactly the mean field value. 



3.2 If =const. limit 

Fixed (/3 means diverging c, for iV — > oo. 



3.2.1 The low temperature case, /3 ^ 1 

For J (3 3> 1, the second term on the right hand side of 
Eq. ( 13 ) vanishes. A lengthy but trivial calculation shows 
, without that term, is larger than zero for 
V < V9 < 1/2, and V < i < 1, where 



< 0. By symmetry, Eq. (13 I 



that Eq. ([13 1 

< 777 < 1, 

we used the fact that y - 
(without the term ~ 1/c) is negative for —1 < m < 0. This 
means there is no phase transition in the thermodynamic 
limit for tp = const, and the system is always in a state of 
maximum magnetization, 777 ~ ±1. 



3.2.2 The high temperature case, /9 <C 1 

For J/3 < 1, we have t = tanh(J/3) 
—771^, and 7/Jnax — "2777, as above. Sett 



Eq. (13 1, we get 



J 13, and j/max 

dlogf/) 



in 



777, = tanh(2777/3J<y5A^) 



(22) 



which means 771 = 1 in the thermodynamic limit. Note, 
that for the case where the coupHng scales as J = Jq/N, 
we get a phase transition, for finite Jo, which is well known 
for the complete graph, i.e. = 1/2. In the context here, 
Lp shifts the phase transition. 

More generally, assume that the coupling scales with 
system size as J = JqN'^, with a < 0, i.e. coupling de- 
creases with size. Assume further that c scales as c = 
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co^f^ with 6 > 0, then ip = cqN^-^ . Using Eq. ([22]) this 
means 

m = tanh(2TO/3JoCoA^"+'') . (23) 

We see that for a — ~b, k = 2/3JoCo can be seen as the 
critical parameter. For a > ~b we get m — ±1, whereas 
for a < —b there is no magnetization, to = 0. 



4 Conclusion 

The crucial observation of this paper is that the summa- 
tion over all topologies in the Ising model on dynamical 
networks is equivalent to re-writing the partition function 
as a sum over all magnetizations. The model - which can 
be seen as a toy model for a variety of socio-economical 
situations - thus drastically reduces in complexity and 
becomes exactly solvable, both for finite size and the two 
possible thermodynamic limits. 

We thank A. AUahverdyan and an anonymous referee 
for various most useful and clarifying comments. Sup- 
ported by Austrian Science Fund FWF Projects P17621 
and P19132. 
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